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ABSTRACT 

Aims. The relation connecting the emitted isotropic energy and the rest-frame peak energy of the vFv spectra of Gamma-Ray Bursts 
(the Amati relation), strictly depends on the cosmological model, so we need a method to obtain an independent calibration of it. 
Methods. Using the Union Supernovae la catalog, we obtain a cosmographic luminosity distance in the y-redshift and demonstrate that 
this parametrization approximates very well the fiducial standard comsomlogical model ACDM. Furthermore, by this cosmographic 
luminosity distance d h it is possible to achieve the Amati relation independent on the cosmological model 

Results. The cosmographic Amati relation that we obtain agrees, in the errors, with other cosmological-independent calibrations 
proposed in the literature. 

Conclusions. This could be considered a good indication in view to obtain standard candles by Gamma-Ray Bursts 
Key words. Gamma rays : bursts - Cosmology : cosmological parameters - Cosmology : distance scale 



1. Introduction 

Supernovae la (SNela) are co nsidered accurate and reliable stan- 
dard candles, (Phil lips et al.h . In recent years, their use as cos- 
mological distance indicators have led to the puzzling discov- 
ery that the Universe is in a phase of accelerated expansion, 
dRiess etal. 1998t iPerlmutter et al. 1999b . This feature has also 
led to the revision of the standard cosmological model, lead- 
ing to wh at is known today as the AC DM concordance model, 
see e.g. (Os triker and Steinhardt 19 95). However it is not pos- 
sible to observe these objects very far in the Universe. The 
most distant Supern ova la was observed at a redshift of z ~ 1.7 
dBenitez et al. 20021) . For this reason, several cosmological anal- 
ysis made by using the various compiled sample of SNela, like 
the Union Catalog, (Kowals ki et al. 20081) . are not able to inves- 
tigate higher redshift regions of the Universe. If we had distance 
indicators at higher redshifts, then we could extend our knowl- 
edge at these unexplored regions. 

One of the possible solutions to this problem could come 
from the Ga mma-Ray Bursts (GRBs) ass umed as cosmologi- 
cal indicators dPiran 20051 : iMeszaros 20061) . GRBs are the most 
powerful explosions in the Universe: this feature allows them to 
be observed at extremely high redshift. The most distant GRB 
observed up to now is at a redshift of ~ 8.1, dTanvir et al. 2009: 
ISalvaterra et al. 20091) . However, GRBs are not standard can- 
dles, since they have no known and well-defined luminosity rela- 
tion. Due to this lack, we have to find another way to use GRBs 
as cosmological beacons. A possible solution could consist in 
finding correlations between photometric and/or spectroscopic 
properties of GRBs themselve s. In the scientif ic literature there 
are several of these relations. (|Schaefer 20061) . One of these is 
the Amati relation, (Amati 2002), which relates the isotropic en- 
ergy emitted by a GRB with the peak energy in the rest-frame 



of the vF(v) electromagnetic spectrum of a GRB. This relation 
has already been widely used to con straining the cosmological 
density parameter (lAmati et al. 20081) . with quite remarkable re- 
sults. However, there is still not a physical link between this cor- 
relation and the mechanisms underlying the production and the 
emission of a GRB. The basic emission process of a GRB is very 
likely not unique, so it is not easy to explain, from a physical 
point of view, such a relation. Recently it has been suggested that 
the Amati relation could depend strongly on the satellite mea- 
surements used fo r detection and the observation of each GRB 
dButler et al. 20071). Howev er this hypothesis has been rejected 
recently, dAmati et al. 20081) . since the relation seems to be ver- 
ified regardless of the satellite considered for the observations 
and detection. 

Although not supported by self-consistent physical moti- 
vations, it is a phenomenological relation which could be ex- 
tremely useful for cosmological considerations. However, a 
problem related with such a relation is that it must be calibrated 
independently of the considered cosmological model. In order to 
compute the energy emitted from an astrophysical object at a cer- 
tain redshift z, we need, as a matter of fact, a measurement of the 
bolometric flux and the distance of the sam e object. For th e first 
quantity, we follow the idea outlined by dSchaefer 2006b : one 
can obtain a very precise measurement of the bolometric fluence 
emitted by a GRB from the observed fluence, the integrated flux 
in the observation time and the spectral model that best fits the 
spectral energy distribution of each GRB. However, the distance 
depends on the considered cosmological model. People usually 
adopt the standard ACDM model, with fixed values of the den- 
sity parameter Q,-. This procedure leads to the so-called circu- 
larity problem when the Amati relation is used to standardize 
GRBs. For this reason we need a cosmology-independent cali- 
bration of the relation. 
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Recently, it was released a calibration with SNela 
data by using different numerical interpolation methods 
dLiang et al. 20081) : the results seem very reliable to address cos- 
mological issues by GRBs. In this work we shall take into ac- 
count a similar analysis: by taking into account SNela data from 
the c osmographic point of view ( for a detailed description see 
e.g. ([Weinberg 19721: IVisser 2 004)). it could be possible to ob- 
tain a calibration of the Amati relation. We will use results ob- 
tained from a cosmographic fit of a sample of SNela extended 
up to very high redshift with the GRBs. The use of the cosmog- 
raphy to deduce the cosmological parameters f rom SNela was 
widely discussed in the literature, (IVisser 200 7a). and the results 
are very close to that attained by other and more accurate analy- 
sis. Recently applications of cosm ographic methods have taken 
into account galaxy cluster s (Capozziello et al. 2004) and GRB, 
dCapozziello & Izzo 20081: IVitagliano et al. 2010h~ but their reli- 
ability drastically fails at high redshifts. Indeed, the estimates 
of the deceleration parameter qo and of the jerk parameter jo 
are usually achieved only at very low redshift and then any ex- 
trapolation could led to shortcomings and misleading results as 
soon as they are extended. However by an appropriate parame- 
terization of the redshift parameter, one can circumvent the prob- 
lem introduc ing a new redshift variable ranging from and 1 
(Visser 2004). Let us consider the following quantity as the new 
redshift variable: 

1 + z 

we obtain that the range of variation is between and 1 . In this 
way, we can derive a luminosity distance by which we can obtain 
the Amati relation suitable for cosmography. 

The layout of the paper is the following: in Sect. 2 we tackle 
the cosmographic analysis considering the SNela Union sample. 
Results will be used to derive the luminosity distance for each 
GRB and then we will fit the cosmographic Amati relation. In 
Sect. 3 a discussion on how to extend the same relation is re- 
ported. We add further 13 GRBs (as of December 2009), com- 
puting the bolometric fluence and the peak energy for each of 
them and after we calculate the cosmographic parameters using 
the new relation. Finally, we calculate the isotropic energy for 
each GRB and then compute the best fit for the considered sam- 
ple of data (Sect. 4). Discussion and conclusions are reported in 
Sec. 5. 



Table 1. SNela cosmographic fit obtained by both the redshift 
variables z and y. 



Parameter 


value z 


■redshift 


Parameter 


value j-redshift 


a 


4242 


± 176 
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4213 ± 216 


b 


0.9593 


t 0.2447 


b 


2.248 ± 0.863 


c 


-0.8201 


± 0.4290 


c 


-0.894 ± 3.576 


d 


0.2722 


± 0.2119 


d 


1.784 ±4.057 
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where dn = c/Hq and Hq, qo, jo and so are known as the 
Hubble constant, the deceleration, the jerk and the snap param- 
eters respectively. In order to obtain accurate measurements of 
the cosmographic parameters, we need to go up to large values 
of the redshift. This goa l can be a chieved by considering large 
data sample as SNela dVisser 2 004) and, eventually, GRBs. 

Here, we are interested in reconstructing the relation d/{z) by 
cosmographic methods in order to test correlations for GRBs . In 
order to achieve this goal, we will use the data sample of SNela 
Union, dKowalski et al. 2 008) consisting of 307 supernovae up 
to redshift z ~ 1.7. By this data sample, it is possible to perform 
a non-linear least-squares fit considering the empirical equation 
given by the distance modulus obtained from the expanded di(z), 
that is: 



H(z) = 25 + — - log \d H [z + -(1 - qo)z 2 
log 10 ( 2 



U, . c 2 k 



-qo- 3ql)z 3 



l r9i r. , ,r 2 , 2c 2 k(l+3qo) 

+ M (2 + 5j0 ~ 2q ° ~ 15q ° + aW 2 



2. Cosmographic analysis 

The main purpose of this work consists in obtaining an Amati 
relation independent of the adopted cosmological model. All we 
need is a formulation of the luminosity distance di as a func- 
tion of the redshift z- These two quantities are linked together 
via the scale factor a(f), which describes the expansion of the 
Universe in a Friedmann-Lemaitre-Robertson- Walker cosmol- 
ogy. This means that we are assuming only homogeneity and 
isotropy but not the specific cosmological model, e.g. ACDM 
model. It is well known that we can obtain the function a(t) 
from the Friedmann equations. These equations can be solved 
only if assumptions are made on dynamics and fluids filling the 
Universe, that is choosing a cosmological model. We will re- 
lax this po ssibility a ssuming only cosmography in the sense de- 
scribed in (We inberg 19721) . Since the evolution of the luminos- 
ity distance is well known for small values of redshift, we can 
consider the power series expansion of the scale factor. This nat- 
urally leads to an expression for the luminosity distanc e in power 
series terms too dVisser 2004 ICapozziello et al. 20081) : 



+ I0joqo-I5ql + S )z 4 + O(z 4 )] 



(3) 



In this work we are not interested in the estimate of the cos- 
mographic parameters but in using cosmography to constrain a 
GRB-energy relation. To this aim, we will use a custom equation 
for the fit of the type 



H(z) = 25 + (5/ log 10) log(az + bz 2 + cz 3 + dz 4 ) , 



(4) 



so we will compute only the parameters a,b,c,d. Once we have 
an estimate of these parameters, we could easily obtain the val- 
ues of the related cosmographic parameters. To obtain a better 
analysis, we can use a robust interpolation method of Levenberg- 
Marquardt type. The results of our data fitting are shown in Table 

m 

T he test of reliability o f the fit has been done with a R 2 - 
test, (IBevington et al. 20021) . whose value is 0.9914. However 
the extension up to high redshift of this function fi(z) shows a 
serious problem: for redshifts greater than ~ 2 the curve grows 
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rapidly, see FigJT] This steep departure is due to the higher- 
order term, i.e. d, which has a decisive influence at high red- 
shift. This fact rules out a priori a possible supernova-calibrated 
fi(z) at high redshift. Such problems can be eliminated if we 
consider a new variable for the re dshift. It has been shown 
dVisser 2007atlVitagliano et al. 20101) that the coordinates trans- 
formation y = z/(l + z) and, consequently, the power series of 
the luminosity distance provides a better extrapolation at high 
redshft, as well as better results for the parameters of the fit. Due 
to this fact, we can perform a cosmographic analysis for the new 
distance modulus /-i(y), in analogy with what has been already 
done for the fi{z). The new expression for the distance modu- 
lus, whic h takes into account the new redshift parameterization, 
becomes (Vitaglian o et al. 20101) : 



5 f 1 
li(y) = 25 + j— jTj j lo S d H + logy - -(qo - 3)y 



1 



c 2 k 



+—(21 - 40o + - n - 5 ) + qo(9q - 2W 

1 c 2 k 
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+2ql - \0ql + s Q ]y 3 + 0(y 4 )] 



(5) 



so we will consider a custom equation for the fit similar to 
the previous one, used for the estimate of the fi(z) parameters. 
The results obtained with a non-linear fit are shown in Table Q] 
while in Figure Q] it is shown the trend of the distance modulus 
for both the redshift variables considered. 



In the following we will consider the formulation for the dis- 
tance modulus in terms of the y-redshift in order to derive a cos- 
mographic Amati relation. 



3. The data sample 

As said in Introduction, in recent years the interest of astrophysi- 
cists and cosmologists has been attracted by the possibility of 
using GRBs as potential distance indicators. This interest is due 
to the fact that most of the GRBs satisfy some correlations be- 
tween photometrical and spectroscopical observable quantities. 
Among th e various exis ting correlations (for a review of these 
see e.g. (ISchaefer 20 06)). the Amati relation seems very attac- 
tive (I Amati 20021) . It relates the cosmological rest-frame v F(v) 
spectrum peak energy E p j with the equivalent isotropic radiated 
energy Ei so . It was discovered based on Bepp oSAX data and the n 
confirmed also for the X-ray flashes (XRFs) (Lamb et al. 2004). 
It seems that it does not work for short GRBs. For this reason 
the relation could be used to discriminate among different GRB 
classes. 

The possible origin of this correlation as due to detector se- 
lection effects seems not consistent, nevertheless the large scat- 
ter in the normalization and the shift toward the Swift detection 
threshold dButler et al. 20071) . A recent study (I Amati et al. 2 009) 
has shown that the different E p j - Ej SO correlations, obtained in- 
dependently from the detectors considered for the observations, 
are fully consistent each other, so the hypothesis of a strumental- 
dependent Amati relation seems to fail. 



Here we are goin g to expand the sample of GRBs reported 
in (lAmati et al. 20091) adding 13 GRBs and obtaining a sample 
consisting of 108 GRBs. Substantially we need to know the red- 
shift z, the observed peak energy E Pt0 i, s of the v F (v) spectrum 
and an estimate of the bolometric fluence S boh for each GRB in 
the sample. To deriv e the bolometri c fluence Sb u we can use the 
method outlined in (ISchaefer 2 006). where from the observed 
fluence and the spectral model, we can obtain an estimation of 
S boi via the following formula: 



>bol 



* obs " 



10 4 /(l+z) 

L/a+z) 



EipdE 



j; i-'Mi- 



(6) 



where 4> is the spectral model considered for the spectral data fit 
and S bs is the fluence observed for each GRB in the respective 
detection band (£„„■„, E max ). In particular, for 6 of the 13 GRBs 
added, we consider a cut-off power-law spectral model w hile for 
the remaining 7 we use a band model dBand et al. 1993b . In the 
Table|2] the spectral data for the 13 GRBs are shown. E p column 
refers to the measured peak energy. To obtain the peak energy in 
the rest frame, we have to take into account the redshift of the 
GRB, then E p j = E p {\ +z). Once we have obtained the estimate 
of S boi for each GRB in the sample, the next step is to estimate 
the isotropic energy from the well-known formula which relates 
the luminosity distance and the fluence, that is 



E iso = 4ndfS b oi(J +z) 



(7) 



Note that the quantity (1 +z) to obtain the value of an observable 
quantity in the rest-frame is equivalent, in the new redshift pa- 
rameterization, to use, instead, the correction 1 /(l— y). The value 
of the luminosity distance which must enter in Eqf7]is what we 
got previously from the cosmographic fit of the SNela. From this 
fit, we obtained an estimate of the function p(y); to go back to 
the luminosity distance, we can use the following formula: 



di(y) = 10" 



(8) 



by which it is possible to compute the value of di(y) for each 
GRB in the sample. 

It is worth noticing that for values of y greater than ~ 2.5, 
the curve fi{y) begins to increase slightly. This fact could lead 
to improper estimates of the isotropic energies emitted by GRBs 
at high redshift. If we consider an analogous curve referred to 
a fiducial standard ACDM cosmological model, we can quanti- 
tatively evaluate this deviation. In Figure [3] it is shown the de- 
viation of the curve fi{y), obtained by the cosmographic fit of 
the SNela and the one obtained by considering a ACDM model 
with values of the density parameters given by Q p = 0.27 and 
Q.A = 0.73. The discrepancy from the fiducial ACDM model 
seems quite small, but it has to be taken into account when we 
will compute the cosmographic Amati relation. 

4. The Cosmographic Amati relation 

At this point we can calculate the parameters of the Amati rela- 
tion for the sample that we constructed previously. This relation 
is a correlation of type Ej SO = aE y p ; ; however if we report it in a 
logarithmic basis, it reduces to the form: 



\og w E iso = A + ylog 10 E p j 



so we can report our sample in a diagram log 10 £, 
and perform a linear fit of the data, with weights given by the 



(9) 

logio E P,i 
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Fig. 1. Trends of the distance modulus for the z-redshift and for the y-redshift. 



Table 2. Data for the 13 GRBs added to the old sample described in dAmati et al. 20091) . Table shows: (1) the name of GRB, (2) the 
spectral model used for the fitting of the spectra, (3) the redshift, (4) the peak energy observed, (5) the softer spectral index, absent 
for the cut-off power law spectral model, (6) the higher spectral index, (7) the observed fiuence and (8) the detector band considered 
for the estimate of the fiuence, (9) the GCN reference for the GRB, where we took the spectral data. 



GRB 


spec model 


z 


E„, (keV) 


a 


B(y) 


S obs (10 6 ergs/cm 2 ) 


band (keV) 


GCN 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


090516 


CPL 


4.109 


190 ± 65 




-1.5+0.3 


15 + 3 


20-1200 


9422 


0907 15B 


CPL 


3.00 


134 ± 56 




-1.1+0.4 


9.3 ± 1.5 


20-2000 


9679 


090812 


CPL 


2.452 


586 ± 243 




-1.03 + 0.07 


26.1 ± 3.4 


15-1400 


9821 


090926B 


CPL 


1.24 


91+2 




-0.13+0.06 


8.7 ± 0.3 


10-1000 


9957 


091018 


CPL 


0.971 


28 ± 16 




-1.53 + 0.59 


1.44 + 0.19 


10-1000 


10045 


091029 


CPL 


2.752 


61.4 ± 17.5 




-1.46 + 0.27 


2.4 + 0.1 


15-150 


10103 


090618 


Band 


0.54 


155.5 ± 11 


-1.26 + 0.06 


-2.50 ± 0.33 


270 + 6 


8-1000 


9535 


090902B 


Band 


1.822 


775 ±11 


-0.696 ± 0.012 


-3.85 + 0.31 


374 + 3 


50-10000 


9866 


090926 


Band 


2.1062 


314 ±4 


-0.75 ± 0.01 


-2.59 ± 0.05 


145+4 


8-1000 


9933 


091003 


Band 


0.8969 


486.2 ± 23.6 


-1.13 + 0.01 


-2.64 ± 0.24 


37.6 ± 0.4 


8-1000 


9983 


091020 


Band 


1.71 


103 ± 68 


-0.93 ± 0.6 


-1.9 + 0.8 


10.4 ± 2.1 


20-2000 


10057 


091127 


Band 


0.49 


36 + 2 


-1.27 + 0.06 


-2.20 ± 0.02 


18.7 ± 0.2 


8-1000 


10204 


091208B 


Band 


1.0633 


124 + 20.1 


-1.44 + 0.07 


-2.32 ± 0.47 


5.8 ± 0.2 


8-1000 


10266 



References: (Sakamoto et al. 2009), (McBreen 2009), (Golenetskii et al. 2009), (Sakamoto et al. 2009b), (Bissaldi & Connaushton 2009), 
dBissaldi 20091) . dBriggs 20090 . dRau 20091). dGolenetskii et al. 2009b|)TdGolenetskii et al. 2009cT) . dBarthelmv et al. 20091) , 

dWilson-Hodge & Preece 20091) . JlvlcBreen 2009bl) 



data errors on both the quantities involved. An /?~-test provides 
an estimation of the reliability of the fit being R 2 = 0.772. This 
is a good value, but not so suitable for our analysis. The results 
of the fit, with errors amounting to a deviation of 3<x, and the 
corresponding covariance matrix are: 

A = 49.154 + 0.306 y = 1.444 ± 0.117 (10) 

0.0136435 -0.005091481 
-0.00509148 0.00197731 J 

A comparison with the results obtained by different inter- 
polation methods (Liang et al. 20081) shows a slight discrepancy 
between the parame ters of the relation . This fact could be due 
to the calibration in dLiang et al. 20081) . It depends on the trend 
traced by SNela, while the cosmographic analysis takes into 
account the corrections due to physical parameters as qo, jo. 
Nevertheless the reason could be another: since the SNela sam- 
ple, use d here to calibrat e the Amati relation, is different from 
that in (Liang et al. 2008]), wher e the authors adopted th e cat- 
alog of 192 SNela discussed in (IWood-vasev et al. 20 07). This 
means that the slight difference in the results could be due to the 
different samples used for the calibration. 



In Figj3] it is shown the plot of the cosmographic Amati re- 
lation. The confidence level curves are calculated as the 3<x de- 
viation from the best fit. Note how the 13 GRBs added to the old 
sample, marked with a circle, are distributed about the best-fit 
curve, indicating that the spectral analysis of these 13 GRBs is 
correct. 



5. Discussion and Conclusions 

The issue to extend the cosmic scale ladder up to medium-high 
redshift is an important questions of modern cosmology. A pos- 
sible way to achieve this goal is to take into account GRBs, the 
most powerful explosions in the Universe. The energy emitted 
by these objects spans about six orders of magnitude. However, 
they cannot be assumed as standard candles in a proper sense. 
Dispite of this lack, the existence of several correlations between 
spectroscopic and photometric observable quantities of GRBs 
allow us to solve in part this problem. The fundamental pre- 
requisite to obtain such relations is to estimate the emitted en- 
ergy in a way independent of the cosmological model. In this pa- 
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Fig. 3. Plot of the cosmographic Amati relation. The line of prediction b ounds represents a deviation of 3<x from the best fit line, the 
thick line. The circle represents the 13 GRBs added to the old sample, (Amat iet al. 20091) . 



55 - 



50 - 



cosmographic fit of the SNela, we have constrained a sample of 
GRBs in a cosmology-independent way so that we have fitted 
a cosmographic Amati relation for GRBs. The results are simi- 
lar to thos e obtained from other analysis performed using other 
methods, dSchaefer 2006t iLiang et al. 20081: 1 Amati 2 002). It is 
important to stress the independence from cosmology and the 
calibration obtained by SNela. In our opinion, this characteristic 
is relevant, from one side, to constrain cosmological models, in 
particular, dark energy models, and, from another side, to check 
the physical validity of the Amati relation. 



35 



0.2 



0.4 



0.6 



o.s 



Fig. 2. Plot of fi(y) computed for a fiducial ACDM cosmolog- 
ical model, the continuous line, and for the reconstructed /j.(y) 
obtained by the cosmographic fit of the SNela, the dashed line. 
Note the slight deviation at very high redshift. 



per, we have considered a relation for the luminosity distance di 
that is independent on the dynamics of the Universe, but, in prin- 
ciple, could work only at small redshift. Although we have use 
a parameterization for the redshift which allows to transform the 
variable z in a new variable y, ranging in a limited interval, we 
have seen that the obtained luminosity distance at high-redshift 
differs slightly from the fiducial model ACDM at high redshift, 
see Fig[3] Nevertheless, since we obtained the curve di(y) by an 
analysis of the SNela Union survey, that extends up to a red- 
shift of ~ 1.7, we achieved an independent estimate at slightly 
higher redshiffl By the way using the di(y) obtained with the 



1 Estimates of the Baryonic Acoustic Oscillations (BAO) performed 
by forthcoming surveys of clusters at intermediate redshift (z ~ 2.5 - 
3.5) will give better approximations of the curve d/(y). 
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